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Abstract 

We present a formulation for potential-density pairs to describe axi- 
symmetric galaxies in the Newtonian limit of scalar-tensor theories of 
gravity. The scalar field is described by a modified Helmholtz equation 
with a source that is coupled to the standard Poisson equation of New- 
tonian gravity. The net gravitational force is given by two contributions: 
the standard Newtonian potential plus a term stemming from massive 
scalar fields. General solutions have been found for axisymmetric systems 
and the multipole expansion of the Yukawa potential is given. In parti- 
cular, we have computed potential-density pairs of galactic disks for an 
exponential profile and their rotation curves. 

Keywords: stellar dynamics - galaxy: kinematics and dynamics - galaxy: 
halo- galaxy: disk - galaxy: structure 

Pacs: 04.50.+h, 04.25.Nx, 98.10,+z, 98.62. Gq, 98.62.Hr 



1 Introduction 

In recent years there has been much interest in the study of dark matter and 
energy in a cosmological and astrophysical context. This has been motived by 
recent independent observational data in the Cosmic Microwave Background 
Radiation (CMBR) at various angular scales [I], type la supernovae 0, and 
the 2dF Galaxy Redshift survey [3], among others. These suggest that O = 
SIa + f^m ~ lj with f^A w 0.7 and Q ro w 0.3, implying the existence of dark 
energy and dark matter, respectively. In this way galaxies are expected to 
possess dark components and, in accordance with the rotation curves of stars 
and gas around the centers of spirals, they might be in the form of halos 0] and 
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must contribute to at least 3 to 10 times the mass of the visible matter The 
origin of these dark components is, however, unknown. 

The theoretical framework to explain the existence of dark components finds 
its origin in theories of elementary particles physics, with the addition of the 
action of gravity. There exist many theories (grand unification schemes, string 
theories, braneworlds, etc) that involve such physics, but scalar-tensor theories 
(STT) of gravity are typically found to represent classical effective descriptions 
of such original theories [Hj- In this way, the scalar fields of these theories are 
the natural candidates to be the quintessence field [7], as a remnant of some 
cosmological function. It has been even suggested that the quintessence field is 
the scalar field that also acts on local planetary scales [8] or on galactic scales 
Moreover, massive scalar fields might account to the dark matter components 
of galaxies in the form of halos. 

Motivated by the above arguments, we have recently studied some STT 
effects in galactic systems ^OJ- We have computed spherical potential-density 
pairs coming from such theories in their Newtonian approximation. In 
the present report we compute general potential-density pairs and take some 
limiting cases. Then, we consider the general axisymmetric case and take, as 
an example, an exponential profile, that typically describes stars on the disk. 
Also, rotation curves of such a model can be directly obtained. 

This paper is organized as follows: in section 2 we present the equations 
of the Newtonian approximation of a general scalar-tensor theory of gravity 
and solutions are given in terms of integrals of Green functions. In section 3, 
we discuss solutions for point-like mass distributions and show how a multipolc 
expansion can be done for the Newtonian limit contribution of the scalar field. 
In section 4, computations are done for general axisymmetric potentials and 
rotation velocities of stars in galaxies. Then, an exponential disk is considered. 
In section 5, we present our conclusions. 

2 Scalar Fields and the Newtonian Approxima- 
tion 

In references ^3 ^] are derived the following differential equations, that are 
valid for a general scalar-tensor theory in its Newtonian limit, 



in which we consider deviations of the scalar field </> from some average value, 
given by the inverse of the Newtonian constant, 1/G; a is a some free constant 
of the theory and m the mass of some boson particle, see details in |1U| . 

The standard Newtonian potential, ip, is obtained when the perturbation is 
set to zero, cf> = and a — 0. Otherwise, the new Newtonian potential is given 



Snap , 



(1) 
(2) 



by 



$N = 1p-^(l>- 



(3) 
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The next step is to find solutions for this new Newtonian potential given a 
density profile, that is, to find the so-called potential-density pairs [TI]. General 
solutions to Eqs. Q and @ can be found in terms of the corresponding Green 
functions 

^ = - fdr s /^-+B.C, (4) 
J |r — r a | 

f o(r )e~ m l r ~ rs 

(6 = 2u dvj-^- : +B.C. , (5) 

J |r — r s | 

and the new Newtonian potential is 

•* = * - - -/"-F^i " " / dr /< t e -tr' + Ba (6) 

The first term of Eq. ®, given by ip, is the contribution of the usual New- 
tonian gravitation (without scalar fields), while information about the scalar 
field is contained in the second term, that is, arising from the influence function 
determined by the Helmholtz Green function, where the coupling a enters as 
part of a source factor. 



3 Point— like masses and multipole expansion 

Now we present solutions for point-like masses and the multipole expansion of 
the potentials which are useful for numerical simulations of galaxies. Substitut- 
ing the following distribution of N point-like masses in Eq. (jSJ 

N 

p(r) = y^m s 6(r-r s ) (7) 

one obtains ^3 HI] 

N N _lr_r \l\ 

with m s a source mass at spatial position r s , and the total gravitational force 
on a particle of mass rrii is 

^F = -V$ A r = m. J a, (9) 

where A = h/mc is the Compton wavelength of the effective mass (m) of some 
elementary particle (boson) determined by specific particle physics models. In 
what follows we will use A instead of mT x . This length can have a range of 
values depending on particular particle physics models. The first right hand 
term in Eq.© is the pure Newtonian part and the second one is the dark 
matter contribution which is of the Yukawa type. There are two limits: On the 
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one hand, for if r ^S> A (or A — * 0) one recovers the Newtonian theory of gravity. 
On the other hand, for if r <C A (or A — > oo) one again obtains the Newtonian 
theory, but now with a rescaled Newtonian constant, G — > Gn(1 + o). There are 
stringent constraints on the possible X—a values determined by measurements 
on local scales [Tl] . 

In the past the above solutions have been used to solve the missing mass 
problem in spirals |151 HB] as an alternative to considering a distribution of dark 
matter. This was done assuming that most of the galactic mass is located in 
the galactic center, and then considering the center to be a point source. In 
our present investigation we do not avoid dark matter, since our model predicts 
that bosonic dark matter produces, through a scalar field associated to it, a 
modification of Newtonian gravity theory. This dark matter is presumably 
clumped in the form of dark halos. Therefore we will consider in what follows 
that a dark halo is axisymmctrically distributed along an observable spiral and 
beyond, having some density profile. Next, we compute the potentials, and 
some astrophysical quantities, for general halo density distributions. 

The standard gravitational potential due to a distribution of mass p(r), 
in a point exterior to the distribution and without considering the boundary 
condition, can be expanded as [17] 

^( r ) = -E E wn qin p+i ' (10) 



21 

1=0 n=-l 



with the expansion 



v~r\ = 4 -E E ^-jfrWW^), (ii) 

where r< is the smaller of |r| and |r s |, and r> is the larger of |r| and |r s |; 
Yi n {6, ip) are the spherical harmonics. 

Here we are interested in the potential outside the distribution of mass, 
consequently, the coefficients of the expansion of ip, and known as multipoles 
are given by 

q ln = ( dr'Y ln (6',^y l p(v'). (12) 



In the case of the scalar field, with the expansion 

7 

exp(— m\r — r s \ 



1 1 tOn=-I 



4 ™E E ii(mr < )k l (mr > )Y&{0',<f/)Y ln (9 1 <p), (13) 



the contribution of the scalar field to the Newtonian gravitational potential can 
be written as 

= E E ^J^inkl(mr)Y ln (9, <p), (14) 
1=0 n—-l 
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where ii{x) and ki(x) are the modified spherical Bessel functions and we have 
defined the multipoles for the scalar field as 



q ln =ml dr f Yi n (e f ,< f /)ii(mr f )p(^). (15) 

4 General axisymmetric density distributions and 
their rotation velocities 

Given that we are interested in axisymmetric galaxies in what follows we will 
only consider the case of this symmetry. Additionally, we use flatness boundary 
conditions (B.C.) at infinity, such that the boundary terms in Eqs. are 
zero. Moreover, regularity conditions must be applied to spatial points where 
the potentials are singular. For axisymmetric systems these conditions mean 
that dtp/dr = d<p/dr = along the symmetry axis. Accordingly, we assume 
that 

PM = j . t > L (16) 

where we define the boundary surface by a single- valued function r = Z(z) with 
z in the interval (zo,zi) and / = V V 2 + z 2 , and the same relation for L with 
(r, z) on the boundary surface. We obtain for tp 

pzi fZ(z) r-oa 

i/)(r,z) = - / dz I r s dr s p(r s ,z s ) / dk 7 (fc?'<)^o(fc? , >) cosfc(z - z s ) , 

(17) 

and for 6 
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/•Zl fZ(z) roo 

^ v r, z) = -a / dz r s dr s p(r s , z s ) / dk I^vr^K^vry) cosk(z-z s ) , 

* Jz Jo Jo 

(18) 

where Io(x) and Kq(x) arc modified Bessel functions of index n — and v 2 — 
k 2 + A -2 . Then, the new Newtonian potential can be obtained using ©. 

An important quantity that characterizes the steady state of astrophysical 
systems is the circular velocity. The circular velocity for a test particle moving 
on the equatorial plane of an axisymmetric galactic system is given by 

v c = r— . (19) 

We may consider as an example that the distribution of mass is 

p{r,z) = Aexp(-r/r )5(z) , (20) 

where tq is the effective radial extension of the system. This distribution of 
mass is typically observed in stars of a disk galaxy. In this case the Newtonian 
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potential is 



$jvO") = ~ A \J r s dr s exp(~r s /r ) J dk I (kr s )K (kr) cos k(z) 

/>OC />OC 

+ / r s dr s exp(-r s /r ) / dk I Q (kr)K Q (kr s ) cos k(z) 

Jr JO 

t>r roc 

+a / r s dr s exp(-r s /r ) / dfc Io(vr s )K (vr) cos k(z) 
Jo Jo 



+a J r s dr s exp(-r s /r ) J dk I (vr)K (vr s ) cos k(z)> . (21) 
Using Eq. (19) the rotation curve is obtained. 

5 Discussion and Conclusions 

We have found potential-density pairs of axisymmetric galactic systems within 
the context of linearized scalar-tensor theories of gravitation. The influence of 
massive scalar fields is given by <f>, determined by Eq. (JTSJ. Circular velocities 
of stars can be directly obtained in the axisymmetric system; by stars we mean 
probe particles that follow the dark halo potential. Specifically, these results 
were used to find potential-density pairs for an exponential profile. In general 
the contribution due to massive scalar fields is non-trivial, see for instance 
Eq. (21), and interestingly, forces on circular orbits of stars depend on the 
parameters A and a in a rather complicated way. This means that even when 
local experiments force a to be a very small number [14| . the amplitude of forces 
exerted on stars is not necessarily very small and may contribute significantly 
to the dynamics of stars. Alternatively, one may interpret the local Newtonian 
constant as given by (1 + a)(</>) -1 , instead of being given by {(f))^ 1 (1 in our 
convention). In this case, the local measurement constraints are automatically 
satisfied, and at scales larger than A one sees a reduction of 1/(1 + a) in the 
Newtonian constant. 

In the past, different authors have used point solutions, Eq. (6), to solve 
the missing mass problem encountered in the rotation velocities of spirals and 
in galaxy cluster dynamics ^3 Ej • These models were used as an alternative 
to avoid dark matter. Indeed, for single galaxies one can adjust the parameters 
(a, A) to solve these problems without the need for dark matter. However, 
these models do not provide a good description of the systematics of galaxy 
rotation curves because they predict the scale A to be independent of the galactic 
luminosity, and this conflicts with observations for different galactic sizes [T%] 
unless one assumes various A's, and hence various fundamental masses, to, one 
for each galaxy size. Such a particle spectrum is not expected from theoretical 
arguments; it represents a considerable fine tuning of masses. This criticism 
would also apply to our models. The purpose of the present investigation was, 
however, not to present a model alternative to dark matter in order to solve 
the missing mass problem, but to compute the influence of scalar-field dark 
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matter distributed in the form of an axisymmetric dark halo. This contribution, 
together with the other components of the galaxy, will give rise to a flat velocity 
curve, in a similar way as in Newtonian mechanics. 
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